Microcanonical phase diagrams of short-range ferromagnets 
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A phase diagram is a graph in parameter space showing the phase boundaries of a many-particle 
system. Commonly, the control parameters are chosen to be those of the (generalized) canonical 
ensemble, such as temperature and magnetic field. However, depending on the physical situation 
of interest, the (generalized) microcanonical ensemble may be more appropriate, with the corre- 
sponding control parameters being energy and magnetization. We show that the phase diagram on 
this parameter space looks remarkably different from the canonical one. The general features of 
such a microcanonical phase diagram are investigated by studying two models of ferromagnets with 
short-range interactions. The physical consequences of the findings are discussed, including possible 
applications to nuclear fragmentation, adatoms on surfaces, and cold atoms in optical lattices. 
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Cooperative effects can lead to remarkable properties 
of many-body systems, and the occurrence of a phase 
transition is a prime example of such an effect. At a phase 
transition, the macroscopic properties of a many-particle 
system change abruptly under variation of a control pa- 
rameter. Typical examples of phase transitions are the 
evaporation of a liquid at temperatures above its boil- 
ing point, or the onset of a spontaneous magnetization 
in a ferromagnet below its Curie temperature. In a ther- 
modynamic description, phase transitions are signaled by 
nonanalyticities of thermodynamic functions like the free 
energy density. For example, when discussing the phases 
and phase transitions of a ferromagnet, the Gibbs free en- 
ergy density g(T, h) as a function of the temperature T 
and the external magnetic field h is considered. The so- 
called phase diagram is obtained by plotting in parameter 
space the nonanalyticities of g, i.e., the points or lines 
in the (T, h) plane at which this function is not infinitely 
many times differentiable. Whenever the parameters T 
and h are varied along a path crossing such a point or 
line of nonanalyticities, the system will undergo a phase 
transition. For a ferromagnet, such a phase diagram has 
a very simple appearance (see Fig. [T]). 
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FIG. 1: Phase diagram of a simple ferromagnet. The bold 
line marks the parameter values in the (T, h) plane at which 
the Gibbs free energy density g{T, h) is nonanalytic. 



This is how a phase transitions appears on the macro- 
scopic, thermodynamic level of description. On a more 
fundamental level, statistical physics provides a micro- 
scopic description underlying the thermodynamic one. 
Under suitable conditions on the interactions, thermody- 
namics can be recovered from the statistical description 
in the thermodynamic limit of infinite number of degrees 
of freedom. In equilibrium statistical physics, statisti- 
cal weights are associated to the various microstates of 
a system, and the choice of these weights depends on 
the physical situation the system is in. Different phys- 
ical situations are described by the so-called statistical 
ensembles. The microcanonical ensemble for example is 
appropriate for the description of an isolated system at 
fixed energy, whereas the canonical ensemble describes a 
system in equilibrium with an infinitely large heat bath 
of temperature T. For a suitable class of short-range in- 
teractions, both ensembles are known to give equivalent 
results in the thermodynamic limit (see |l| for details). 
For long-range interactions, however, equivalence may be 
violated, in particular whenever a discontinuous phase 
transition takes place in the canonical ensemble 0]. 

But also in the case of short-range interactions when 
different ensembles yield equivalent results in the thermo- 
dynamic limit, an important difference remains, which 
we would like to discuss, taking again a simple ferro- 
magnet as an example: Canonically, the temperature T 
and the external magnetic field h are the relevant con- 
trol parameters. Microcanonically, however, this is no 
longer true. In a (generalized) microcanonical ensemble, 
the energy (density) e and the magnetization (density) 
m are the natural control parameters corresponding to 
T and h [13] ■ Working in this microcanonical ensemble, 
one would naturally ask: What does the phase diagram 
corresponding to the one in Fig. [1] look like in the (e, m) 
plane? This diagram then could readily answer the ques- 
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tion whether, upon variation of e and m along a certain 
path, the system undergoes a phase transition or not. 
Remarkably, such a microcanonical phase diagram looks 
very different from its canonical counterpart. Moreover, 
although the microcanonical ensemble is the most fun- 
damental one among the statistical ensembles, little is 
known about these issues. To some extend this is also 
due to the fact that studies in the microcanonical ensem- 
ble, either analytically or numerically, are typically more 
demanding than in the canonical ensemble. 

In this Letter, we present two case studies which serve 
to illustrate the general properties of microcanonical 
phase diagrams of short-range interacting ferromagnets. 
The first example is the spherical model with nearest- 
neighbor interactions on a d-dimensional hypercubic lat- 
tice, for which analytical results are presented. Sec- 
ond, numerical results are reported for the Ising model 
with nearest-neighbor interactions on a two-dimensional 
square lattice. Remarkably, both models show two dis- 
tinct transition lines in the (e, m) plane which, for fixed 
magnetization m, may be crossed upon variation of the 
energy e. In the conclusions, these unexpected results 
are discussed, in particular as what regards physical re- 
alizations of an ensemble where both e and m are fixed. 

Spherical model. — Introduced by Berlin and Kac Q in 
1952, this model was constructed to show a ferromag- 
netic phase transition while being exactly solvable. The 
degrees of freedom (7^ S R are associated to the sites 
of a d-dimensional hypercubic lattice. The energy of a 
microstatc a = (di, . . . , aj^) of N degrees of freedom is 

- -J^a.CTj -/i^a„ (1) 
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where J > is a coupling constant determining the 
strength of the exchange interaction. The angular brack- 
ets denote a summation over all pairs of nearest neighbors 
on the lattice. In addition, the ai are required to satisfy 
the spherical constraint J2i'^i — which accounts for 
the model's name. In the canonical ensemble, the spheri- 
cal model is exactly solvable in the thermodynamic limit 
for arbitrary spatial dimension d, and a transition from 
a ferromagnetic phase at low temperatures to a param- 
agnetic phase at high temperatures occurs for d^ i. 

Starting point for a calculation in the microcanonical 
ensemble is the density of states as a function of energy 
e and magnetization m. 



nN{£,m)= [ daS[Ne-H{a)] 

x6[Nm-M{a)]sfN-J2<^^), (2) 



where AI{a) — yields the total magnetization of 

a microstate. An analytic calculation of 57 at is reported in 
[4j|, but the saddle point analysis proposed in that paper 
works only within a certain range of e and m values. This 
can be seen by starting from Eqs. (25) and (26) of Ref. 
^ and, by explicitely performing two of the integrations, 
rewriting the density of states for large N as 
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A derivation of this result will be given elsewhere. Due 
to the logarithm in ([3|) , the integrand of the z-integration 
has two branch cuts on the real line. Apart from these 
branch cuts, the integrand is holomorphic, and the con- 
tour of integration can be deformed freely, as long as it 
does not cross the cuts. For an asymptotic evaluation of 
the integral in the large- iV-limit by means of the method 
of steepest descent, the path of integration in the complex 
plane is deformed such that its imaginary part becomes 
zero. In this limit, the value of the integral is given by 
the integrand of the z-integration in ([3]) evaluated at the 
maximum along that path. Depending on the values of 
e and m, this maximum may either be a saddle point of 
the exponent in ^ (and in this case the analysis in [4] is 
valid), or located at one of the end points of the branch 



cuts. The transition between these two types of behav- 
ior accounts for nonanalyticities of the microcanonical 
entropy in the thermodynamic limit [l8| . 



s(e,m)= lim — In riAr(£,m) 



(4) 



and hence for the occurrence of phase transitions [If 
Note that, contrary to what has been conjectured in 
0|, the entropy s of the spherical model is found to be 
a concave function on its entire domain. By means of 
an asymptotic analysis of Eqs. ([3]) and (j4|), the values of 
e and m for which s(£, m) becomes nonanalytic can be 
computed, yielding 
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with 

Plotting the two curves e±(m) in the (e,™) plane, the 

microcanonical phase diagram of the spherical model is 

obtained (see Fig. [2] for a plot of the d = 3 case) . 
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FIG. 2: (Color online) Microcanonical phase diagram of the 
spherical model on a three-dimensional cubic lattice. The mi- 
crocanonical entropy s is defined only within the gray shaded 
region in the (e, m) plane. Within each of the three gray 
shaded subregions, the entropy is analytic, but not so on their 
boundaries which are given by Eqs. (O and 

This microcanonical phase diagram looks remarkably 
different from its canonical counterpart in Fig. [1] and 
even an experienced statistical physicist, we suspect, 
would have had problems predicting its shape. Varying, 
for example, the energy e while keeping fixed the magne- 
tization at any value of m, one typically crosses two tran- 
sition lines in the phase diagram, therefore observing two 
phase transitions, signaled by kinks in the specific heat. 
Similarly, four transition lines are crossed upon variation 
of the magnetization while keeping the energy e fixed at 
any value — d < £ < e+(0), while two transition lines are 
crossed for energies e+(0) < £ < e_(0). 

The dashed line in Fig. [2] corresponds to a transition 
from a ferromagnetic to a paramagnetic phase, and the 
region to the left of this line is the coexistence region. 
When switching to the canonical ensemble by means of 
a Legendre-Fenchel transform, the entire coexistence re- 
gion is mapped onto the transition line in the canonical 
phase diagram. The solid line in Fig. [21 in contrast, has 
no counterpart in Fig.[Tl mainly due to the fact that this 
transition occurs at negative microcanonical inverse tem- 
peratures ds/de. So what kind of phase is then found to 
the right of this second transition line? At least for van- 
ishing magnetization m one can argue that, upon cross- 
ing this lineja transition to an antiferromagnetic phase 
takes place [20|. For m 7^ 0, however, an interpretation 
of the transition is more difficult since m is not an order 
parameter of the antiferromagnetic transition. 

Two-dimensional Ising model. — The Ising model is ar- 
guably the most studied model in the theory of phase 
transitions, serving as a test case also for our aim of 
computing the microcanonical phase diagram. Its energy 
function is formally equivalent to that of the spherical 



model on a two-dimensional square lattice, 

i/(tTl, . . . ,(Tjv) = O-^fTj - /l^CTj, (7) 

but the degrees of freedom a-i G {—1, +1} take on only 
discrete values. Again, the angular brackets denote a 
summation over all pairs of nearest neighbors on the lat- 
tice. For vanishing external field h and in the thermody- 
namic limit, this model is known to undergo a phase tran- 
sition from a ferromagnetic phase at low temperatures to 
a paramagnetic phase at high temperatures, taking place 
at a critical inverse temperature /3c = ln(l -I- ■\/2)/2. The 
analytic solution for the free energy density was obtained 
by Onsager in 1944 Q. For h ^ 0, however, no analytic 
solution is known. In the microcanonical framework, this 
corresponds to the fact that analytic results exist only for 
some regions in the (£, m) plane, but not for all. Conse- 
quently, we will resort to numerical methods in order to 
compute the microcanonical phase diagram. 

We use a Monte Carlo histogram method discussed in 
in the context of the Ising model with fixed magne- 
tization. For square Ising systems composed of N = 
lattice sites, we compute the density of states r^Af (£, wo) 
for a clamped value of the magnetization m = mo- In a 
microcanonical analysis, quantities of interest are directly 
derived from f2jv(£, mg) or, equivalently, from the mi- 
crocanonical entropy SAr(£,mo) = ln07v(£, toq)/-^- The 
microcanonical specific heat at fixed magnetization m = 
TOO) for instance, is given by c{e) — —[ds/deY/[d'^s/de'^]. 
It is this quantity that we use for our investigation of 
the microcanonical phase diagram of the two-dimensional 
Ising model, mainly by looking for peaks in the specific 
heat which can be viewed as finite-system precursors of 
nonanalyticities occurring in the thermodynamic limit. 

Fig. [3] summarizes our results for systems composed of 
100 X 100 spins. When looking at the microcanonical spe- 
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FIG. 3: (Color online) (a) Microcanonical specific heat as 
a function of energy e for fixed values of the magnetization 
m = mo. From top to bottom: mo = 0,0.05,0.1,0.125,0.15. 
The peaks reveal the presence of a phase transition line at 
relatively high energies. The data shown have been obtained 
for systems composed of 100 x 100 spins, (b) Resulting micro- 
canonical phase diagram of the two-dimensional Ising model. 
The dashed line is the line of spontaneous magnetization sep- 
arating the ferromagnetic phase at low energies from the para- 
magnetic phase. The solid line is new and should be compared 
with the corresponding line in the microcanonical phase dia- 
gram of the spherical model in Fig. (2] 
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cific heat at a certain (for practical reasons not too large) 
value of the fixed magnetization, we observe two peaks. 
One peak occurs at lower energies, signaling the crossing 
of the coexistence line which separates the ferromagnetic 
phase from the paramagnetic phase (not shown) 
A second peak is observed for much larger energies, as 
shown in Fig. [3{a). The positions of these peaks shift 
to smaller energies when the fixed value of the magneti- 
zation increases. We have plotted the peak positions in 
the (e, m) plane in Fig.[3jb), yielding the microcanonical 
phase diagram of the Ising model. The phase diagram 
of the Ising model strikingly resembles the correspond- 
ing phase diagram of the spherical model shown in Fig. 
[2j In both cases there is a range of magnetizations for 
which lines at fixed energy cross two different transition 
lines. Note that for the Ising model, we are not able to 
determine whether the new line extends all the way down 
to the ground state, as this line rapidly closes in on the 
boundary of the microcanonical entropy's support. 

In order to check the robustness of our findings for the 
Ising model, we also studied larger systems with up to 
600 X 600 spins. We observe that the peaks increase in 
height and get sharper for increasing system sizes, as ex- 
pected for a phase transition. The positions of the peaks 
shift slightly towards larger values of e when increasing 
the system size, but this shift is so small that it remains 
within the thickness of the line on the scale of Fig. [3l[b). 
For very large systems, it is extremely difficult to obtain 
the high quality data needed for a microcanonical anal- 
ysis, and we were not able to make a quantitative study 
of the change in peak height and position. 

Discussion. — The phase diagrams discussed above rep- 
resent a physical situation in which both energy e and 
magnetization m can be controlled externally. Consider- 
ing the spherical model or the Ising model properly as 
models of ferromagnets, control of the energy may well 
be imagined in an experimental set-up energetically iso- 
lated from the environment. Direct control of the magne- 
tization, however, appears difficult — if not impossible — 
to achieve. But ferromagnetic spin models have a wide 
range of applications, going well beyond the modeling 
of ferromagnetic materials, both in classical and quan- 
tum physics. Following Lee and Yang, the Ising model 
can be mapped onto a lattice gas, in which the magne- 
tization within the first model formally corresponds to 
the particle density within the latter Control of 



the particle density is of course an experimentally real- 
istic scenario, and in this situation the microcanonical 
phase diagrams of ferromagnetic models can provide rel- 
evant information. Examples include the Ising model as a 
model of nuclear matter fragmentation or as a model 
of adatoms on a crystal surface [l^. Cold atoms in an 
optical lattice are another possible experimental realiza- 
tion: After switching off the cooling, total energy and 
number of atoms are conserved to a very good degree, 
rendering appropriate a description as a lattice gas in the 



microcanonical ensemble. The interactions between the 
atoms can be tuned via Feshbach resonances, allowing to 
realize, among others, Ising-type interactions [3]. 

A comment is in order on the short-range nature of the 
interactions in the two examples discussed. Although we 
believe that the qualitative behavior of the microcanon- 
ical phase diagram should not be restricted to nearest- 
neighbor interactions, it surely does not extend to fer- 
romagnets with long-range interactions. This becomes 
obvious when considering for example the mean-field </3^ 
model. Although this model undergoes a ferromagnetic 
transition, a calculation of the microcanonical entropy 
s(e, m) as a function of energy and magnetization yields 
a smooth function [l^, [l^ . Therefore, the microcanon- 
ical phase diagram of the mean-field tp"^ model in the 
(£,m) plane does not show any transition lines at all. 
This can be seen as a consequence of the nonequivalence 
of the microcanonical and the canonical ensemble in this 
long-range interacting model. Apart from short-range in- 
teractions, we also expect an upper bound on the energy 
per particle to be essential for the observed behavior. 
This is usually the case for spin models, often allowing 
for negative microcanonical temperatures to occur. 

Conclusions. — We have computed microcanonical 
phase diagrams in the parameter space of energy and 
magnetization for two ferromagnetic models. The dia- 
grams look remarkably different from the corresponding 
canonical ones in the (T, h) plane, with the consequence 
that, when controlling e and m in a microcanonical set- 
ting, the physical behavior differs significantly from the 
canonical situation in which T and h are controlled. For 
both models investigated, the microcanonical phase di- 
agrams are qualitatively similar to each other, and we 
expect this to extend to short-range ferromagnetic spin 
models more generally. Finally, we have pointed out 
physical applications of microcanonical phase diagrams 
within the lattice gas interpretation of the Ising model, 
including nuclear fragmentation, adatoms on surfaces, 
and cold atoms in optical lattices. 
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